clear all

import excel "../../data/raw data/Hiroshima_Shisei_Yoran_Sengo/buildingbydistance_1946_changerate.xlsx",  firstrow

* 全壊率＝Total destruction rate
* 変化率45NOV = Employment change rate before and after the bombing
gen 全壊率2=全壊率^2
reg 変化率45NOV 全壊率 全壊率2

#d;
twoway (scatter 変化率45NOV 全壊率, sort color(black) msize(small)) 
(qfit 変化率45NOV 全壊率, sort color(cranberry) lwidth(medthick) range(0.0 100.0)),
	scheme(s1color) plotregion(lwidth(none) ilwidth(none))
   legend(off)
	xtitle("Percentage of totally destroyed buildings")
	ytitle("Establishment change rate due to the bombing") ylabel(, angle(horizontal));
	graph export "../../output/figure/Estchange_bomb_tdest.pdf", as(pdf)   replace ;
#d cr